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The Belle Collaboration has recently reported a study of the decays B — >■ Ds\{2b2i6)'^ D^*' and 
has given also estimates of relevant ratios between branching fractions of decays B — > D^*' D*j 
providing important information to check the structure of the _Dso(2317), Dsi(2460) and _Dsi(2536) 
mesons. The disagreement between experimental data and Heavy Quark Symmetry has been used 
as an indication that _D3o(2317) and -Dsi(2460) mesons could have a more complex structure than 
the canonical cs one. We analyze these ratios within the framework of a constituent quark model, 
which allows us to incorporate the effects given by finite c-quark mass corrections. Our findings 
are that while the -Dsi(2460) meson could have a sizable non-qq component, the Z)*o(2317) and 
Dsi(2536) mesons seem to be well described by a pure qq structure. 
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I. INTRODUCTION 
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i?-factoricS; through e+e~ collisions, tuned to different 
bottomoniuni energy resonances, have become a source 
of data on heavy hadrons. Bottonionium states decay 
most of the cases to BB pairs, and these B mesons decay 
subsequently into charmed and charmless hadrons via the 
weak interaction. 

Nonleptonic decays of B mesons, described at quark 
level by an effective four-quark interaction h — > ccs, have 
been used to search for new charmonium and charmed- 
strange mesons and to study their properties in detail. 
Within the charmed-strange sector, the BaBar Collabo- 
ration found, in the inclusive I?j~7r° invariant mass dis- 
tribution from e^e~ annihilation data, the narrow state 
£'*o(2317) [1]. The CLEO Collaboration, aiming to con- 
firm the previous state, observed its doublet partner 
Dsi(2460) in the D*+t:° final state [2]. However, the 
properties of these states were not well known until the 
Belle Collaboration observed the B ^ DD*f^{2317) and 
B -^ DD si{2460) decays [3]. 

First observations of the B -^ Z)(*'Dsi(2536) de- 
cay modes have been reported by the BaBar Collabo- 
ration [4, 5] and an upper limit on the decay B'^ — >■ 
D*~Dsi{2536)~^ was also obtained by the Belle Collabo- 
ration [6] . The most recent analysis of the production of 
Z?si(2536)"'" in double charmed B meson decays has been 
reported by the Belle Collaboration in Ref. [7]. Using 

the latest measurements of the B — > D''*'^D*j branching 
fractions [8] they calculate the ratios 
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= 0.10 ±0.03, 
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B {B ^ DD, i{2A60)) 
B{B -^ DD*) 

B{B ^ D*D, i{2m0)) 
B{B -^ D*D*) 



0.44 ±0.11, 



0.58 ±0.12. 
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In addition, the same ratios are calculated for B — > 
D^*''D si {2536)~^ decays, using combined results by the 
BaBar [5] and Belle [7] Collaborations 
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BiB ^ DDsi{2536)) 
B{B -^ DD*) 

B{B^ D*Dsi{2536)) 
B{B^ D*D*) 



= 0.049 ±0.010, 



= 0.044 ±0.010. 



(2) 



/^From a theoretical point of view, this kind of de- 
cays can be described using the factorization approxi- 
mation [9] . This amounts to evaluate the matrix element 

which describes the B — >■ D^*^D]*j weak decay process 
as a product of two matrix elements, one for the B weak 
transition into the D^*") meson and the second one for the 
weak creation of the cs pair which makes the D*j meson. 
The latter matrix element is proportional to the corre- 
sponding D),j meson decay constant and therefore these 
processes can give information about the Dgj structure. 
The 15*0(2317) and L',si(2460) states belong to the 
doublet j^ = i predicted by Heavy Quark Symme- 
try (HQS). They have surprisingly light masses and are 
located below DK and D*K thresholds, respectively. 
This implies that these states are narrow. Neither quark 
models nor Lattice QCD calculations [10] can give a 
satisfactory explanation of this doublet in the charmed 
and charmed-strange sectors. This fact has stimulated 
a fruitful line of research, suggesting that their struc- 
ture is much richer than what one might guess assum- 
ing the qq picture, which has proved quite successful in 



other sectors. Some authors have suggested a qqqq struc- 
ture [11, 12] or mixing between the usual qq structure 
and four quark [13]. Also there have been suggestions 
that these states might be molecular states. However, 
one must mention that the B{B —>■ D*q{2317)K) is the 
same order of magnitude as B{B -^ DgK) and at least a 
factor of 2 larger than the B{B -> D,si{2460)K) branch- 
ing fraction [14], in contrast with the naive expectation 
that decays within the same spin doublet would have sim- 
ilar rates. This fact suggests that the £'*o(2317) meson 
may be a conventional 0+ cs state and the I?si(2460) has 
a more complex structure. 

In any case the study of nonleptonic decays can help 
to shed light on the structure of the P-wave open charm 
mesons 

The D).j meson decay constants are not known experi- 
mentally except for the ground state, Dg, which has been 
measured by different collaborations. Another way to 
study DIJ mesons that does not rely on the knowledge of 



M 



their decay constants is through the decays B^ — ?> D),j' 
where M is a meson with a well known decay con- 
stant. However, the experimental study of these pro- 
cesses is currently difficult leaving, for the time being, 
the B — > D'^*'>D^*j decay processes as our best option to 
study D),j meson properties. 

According to Ref. [15], within the factorization approx- 
imation and in the heavy quark limit, the ratios Rdo and 
Rdi can be written as 



Rdq — Bd'o — 
Rdi = Rd'I = 



/d*o(2317) 
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/Dsi(2460) 
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expect they could also play an important role in this case. 

We work within the framework of the constituent 
quark model described in Ref. [23] and which is widely 
used in hadronic spectroscopy [23, 24]. It has re- 
cently been applied to mesons containing heavy quarks 
in Refs. [25, 26]. It has also been used successfully in 
the description of scmilcptonic B — ^ D** and Bs ^- -D** 
decays in Ref. [27], which provides us with confidence 
that the model describes well the weak decay transitions 
B -^ D^*^ which is one of the terms in the calculation of 
B -^ D^*^d[*J decays. 

The paper is organized as follows: In Sec. II, we intro- 
duce the constituent quark model and discuss its pre- 
dictions in the charmed and charmed-strange sectors. 
Sec. HI is devoted to explain the theoretical framework 
through which we calculate the B — )■ D'^*^D*j decays. 
Finally, we present our results in Sec. IV and give the 
conclusions in Sec. V. 



II. CONSTITUENT QUARK MODEL 

Constituent light quark masses and Goldstone-boson 
exchanges coming from the spontaneous chiral symme- 
try breaking of the QCD Lagrangian, together with the 
perturbative one-gluon exchange (OGE) and the nonper- 
turbative confining interaction are the main pieces of po- 
tential models. Using this idea, Vijande et al. [23] de- 
veloped a model of the quark-quark interaction which is 
able to describe meson phenomenology from the light to 
the heavy quark sector. 

A consistent description of light, strange and heavy 
mesons requires an effective scale-dependent strong cou- 
pling constant. We use the frozen coupling constant of 
Ref. [23] 



where the phase space effects arc neglected because they 
are subleading in the heavy quark expansion. Now, in the 
heavy quark limit one has fo'^ = Jd^i and fo. = Id* 
and so one would predict Rdo ~ Rdi- Moreover, there 
are several estimates of the decay constants, always in the 
heavy quark limit [16-18], that predict for P-wave jq = 
1/2 states similar decay constants as for the ground state 
mesons (i.e. Dg and -D*), and very small decay constants 
for P-wave, jq = 3/2 states. Then, these approximations 
lead to ratios of order one for D*o(2317) and Dsi{2460), 
in strong disagreement with experiment. 

Leaving aside that the factorization approximation has 
been recently analyzed in Refs. [19-21] finding that it 
works well in these kind of processes, we will concen- 
trate in the influence of the effect of the flnite c-quark 
mass in the theoretical predictions. As found in Ref. [22], 
l/jTiQ contributions could give large corrections to var- 
ious quantities describing B — >■ D** transitions and we 



as{p,j) 



ao 



In 






(4) 



where fiij is the reduced mass of the qq pair, and ckq, 
fiQ and Aq are parameters of the model determined by a 
global fit to all meson spectrum. 

In the heavy quark sector chiral symmetry is explicitly 
broken and Goldstone-boson exchanges do not appear. 
Thus, OGE and confinement are the only interactions 
remaining. 

The one-gluon exchange potential is generated from 
the vertex Lagrangian 



Cqqg = i^/4mXs^pJf,G'^X''1p, 



(5) 



where A^ are the 5C/(3) color matrices and G^ is the 
gluon field. The resulting potential contains central, ten- 
sor and spin-orbit contributions given by 
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tors. Note the contact term of the central part of the 
one-gluon exchange potential has been regularized as fol- 
lows 
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Although there is no analytical proof, it is a general 
belief that confinement emerges from the force between 
the gluon color charges. When two quarks are separated, 
due to the non-Abclian character of the theory, the gluon 
fields self-interact forming color strings which bring the 
quarks together. 

In a pure gluon gauge theory the potential energy of 
the qq pair grows linearly with the quark- antiquark dis- 
tance. However, in full QCD the presence of sea quarks 
may soften the linear potential, due to the screening of 
the color charges, and eventually leads to the breaking 
of the string. We incorporate this feature in our confine- 
ment potential as 
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where as controls the mixture between the scalar and 
vector Lorentz structures of the confinement. At short 
distances this potential presents a linear behavior with 
an effective confinement strength, a = —Ucfic (AJ • A^), 
while it becomes constant at large distances. This type 
of potential shows a threshold defined by 



l^thr = {-flc + A}(A? • A^). 



(9) 



No qq bound states can be found for energies higher 
than this threshold. The system suffers a transition 
from a color string configuration between two static color 
sources into a pair of static mesons due to the breaking 
of the color fiux-tube and the most favored subsequent 
decay into hadrons. 



Quark masses 


m„ (MeV) 


313 




rus (MeV) 


555 




mc (MeV) 


1763 




mt (MeV) 


5110 


OGE 


Qo 


2.118 




Ao (fm-i) 


0.113 




A*o (MeV) 


36.976 




fo (fm) 


0.181 




fg (fm) 


0.259 


Confinement 


Oc (MeV) 


507.4 




Mc (fm-^) 


0.576 




A (MeV) 


184.432 




fls 


0.81 



TABLE I. Quark model parameters. 



Among the different methods to solve the Schrodinger 
equation in order to find the quark-antiquark bound 
states, we use the Gaussian Expansion Method [28] which 
provides enough accuracy and it simplifies the subsequent 
evaluation of the decay amplitude matrix elements. 

This procedure provides the radial wave function solu- 
tion of the Schrodinger equation as an expansion in terms 
of basis functions 



Rcir)^ E 



€iir). 
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where a refers to the channel quantum numbers. The 
coefficients, c", and the eigenvalue, E, are determined 
from the Rayleigh-Ritz variational principle 
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(T«„-£;iv«jc^ + E K 



(11) 



where T'",^, iV",^ and V^" are the matrix elements of 
the kinetic energy, the normalization and the potential, 
respectively. T",„ and N",^ are diagonal, whereas the 
mixing between different channels is given by V^,'^ . 

Following Ref . [28] , we employ Gaussian trial functions 
with ranges in geometric progression. This enables the 
optimization of ranges employing a small number of free 
parameters. Moreover, the geometric progression is dense 
at short distances, so that it enables the description of 
the dynamics mediated by short range potentials. The 



fast damping of the gaussian tail does not represent an 
issue, since we can choose the maximal range much longer 
than the hadronic size. 

Model parameters fitted over the whole meson spec- 
tra [23, 25] are shown in Table I. 

The model described above is not able to reproduce the 
spectrum of the P-wave charmed and charmed-strange 
mesons. The inconsistency with experiment is mainly 
due to the fact that the mass splittings between the 
i:>;o(2317), D^i(2460) and D,i{2536) mesons are not weU 
reproduced. The same problem appears in Lattice QCD 
calculations or other quark models [10]. 

In order to improve these mass splittings we follow the 
I 



proposal of Rcf. [29] and include one-loop corrections to 
the OGE potential as derived by Gupta et al. [30]. This 
corrections shows a spin-dependent term which affects 
only mesons with different flavor quarks. 

The net result is a quark- ant iquark interaction that 
can be written as: 

V{n,) = VoG^in,) + Vco^{n,) + ^oge°'(^^-), (12) 

where Voge and Vcon were defined before and are 
treated non-perturbatively. V'qq^"^ is the one-loop cor- 
rection to OGE potential which is treated perturbatively. 
As in the case of Voge and Vcon, ^oge°'' contains cen- 
tral, tensor and spin-orbit contributions, given by [29] 
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where Cf = 4/3, Ca = 3, 60 = 9, 7s = 0.5772 and the 
scale ^ ^ 1 GeV. 

Table II shows the masses of well established charmed 
and charmed-strange mesons predicted by the con- 
stituent quark model and those including one-loop cor- 
rections to the OGE potential. 

The masses predicted for the 0~ and 1^ states agree 
with the experimental measurements in both sectors. 
These states can be identified with the members of the 
lowest lying j^ = ^ doublet predicted by HQS. 

The doublet jq — ^ , which corresponds to the 2+ 
state and one of the low lying 1+ states, is in reasonable 
agreement with experiment. 

The charmed and charmed-strange 0+ states are sen- 
sitive to the one-loop corrections of the OGE potential 
which bring their masses closer to experiment. However, 
the spin dependent corrections to the OGE potential, are 



I 

not enough to solve the puzzle in the 1+ sector. A pos- 
sible explanation for the low mass of this state has been 
proposed in Ref. [26] were the 1+ cs states are coupled to 
a csnn tctraquark structure. One finds that the J^ = 1"*" 
1581(2460) has an important non-qq contribution whereas 
the J^ = 1+ I?si(2536) is almost a pure qq state. 



III. NONLEPTONIC B -^ D'-'^D]*/ DECAYS 



We give account of the nonleptonic decays B — > 
D^*'^D^*] with £)(*) the D or D* mesons, and d[*] the 
mesons D^o (23 17), £>si(2460) and D,i{2hm). These de- 
cay modes involve a b ^- c transition at the quark level. 



This work (os) 

This work (q^) 

Exp. 



Charmed 



if = 1/2' 



if - 1/2^ 



= 3/2+ 



0" 

1896 

1896 

1867.7 ±0.3 



1" 

2017 

2014 

2010.25 ±0.14 



0+ 

2516 

2362 

2403 ± 38 



1+ 

2596 

2535 

2427 ± 36 



1+ 

2466 

2499 

2423.4 ±3.1 



2+ 

2513 

2544 

2460.1 ±4.4 



if = 1/2- 



Charmed-strange 
if - 1/2+ 



j^ = 3/2^ 



This work (os) 

This work (q^) 

Exp. 



0" 1" 0+ 1+ 1+ 2+ 

1984 2110 2510 2593 2554 2591 

1984 2104 2383 2570 2560 2609 

1969.0 ±1.4 2112.3 ±0.5 2318.0 ±1.0 2459.6 ± 0.9 2535.12 ± 0.25 2572.6 ± 0.9 



TABLE II. Masses of well established charmed and charmed-strange mesons predicted by the constituent quark model {cts 
and those including one- loop corrections to the OGE potential (a^). 



governed by the effective Hamiltonian [9, 31, 32] 



H, 



cff. 



^^ {Vcb [Ci(m') Qf ± C^ifi') Qf ] + h.c.} , 



V2 



(14) 

in which penguin operators have been neglected. In 
Eq. (14), Ci{fj,') and C2(/i') are scale-dependent Wilson 
coefficients, being //' ~ nib the appropriate energy scale 
in this case. The Qf' and Qj^ are local four-quark oper- 
ators given by 

Qf = K: [V^c(0)7m(I- 75)^-6(0)] [^.(0)7^(X-75)V'c(0)] , 

(15) 



where Gp = 1.16637(1) x 10"^ GcV"^ is the Fermi decay 
constant [8], A(a,6,c) = (a + & — c)^ — 4a6, Vcb and Vcs are 
the corresponding Cabbibo - Kobayashi- Maskawa matrix 
elements, for which we take Vbc — 0.0413 and Vcs = 
0.974. Pb, Pom and P^(*) are the meson momenta. 

T-LaisiPBiPoi')) is the hadron tensor for the B~ -^ D(*)° 
transition and i-L"'^ [P {,)) is the hadron tensor for the 

vacuum — > D*j transition. We have 



for a scalar or pseudoscalar meson, and 



(18) 
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We show in Fig. 1 the schematic representation of the 
B^ -^ D''*^^D*j decay given by the two local four- 
quark operators, diagram di for Qf' and diagram ^2 for 
Qf'. Factorization approximation is implicit in diagram 
di , which amounts to evaluate the hadron matrix element 
of the effective Hamiltonian as a product of quark-current 
matrix elements. Fierz reordering of diagram ^2 leads to 
the same contribution but for a colour factor. The full 
transition amplitude can be evaluated with the Qf part 
of the Hamiltonian but with an effective coupling given 
by 



with Nc = 3 the number of colors. 
The decay width is given by 



(16) 
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(17) 



for a vector or an axial- vector meson. In previous equa- 
tions, fpi') is the D),/ meson decay constant. 

As in Refs. [27, 31], the contraction of hadron ten- 
sors, 'Hc«/3(Pb^, -?£!(.) )'H"^(-Pp(.)), can be easily written 

in terms of helicity amplitudes for the B~ — > D^*^ tran- 
sition [9], so that the decay width is given as 
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for D*j a pseudoscalar or scalar meson, and 
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(21) 
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FIG. 1. Schematic representation of B decay into D'^'^'^D^j 



for Dlj a vector or axial- vector meson. The different 
Hrr are defined in Ref. [31] and evaluated in this case 
at q^ = m? (,, . The expressions of the decay constants 

needed are given in the Appendix. 



IV. RESULTS 

In this section we present our results for the ratios 
in Eqs. (1) and (2). We have calculated them using 
the factorization approximation. As we have already 
mentioned, the factorization approximation and heavy 
quark limit predicts ratios of order one for _D*o(2317) and 
D,i(2460) mesons, and negligible for D,si(2536) [15, 33]. 

^From Eqs. (20) and (21) we arrive at 
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and the same for Rd'o-, Rd*i and Rd'V but replacing 
the meson D by the meson D* . 

Using experimental masses we obtain the ratios 



Ron = 0.9008 x 
Rd'q = 0.7166 X 



/n*o(2317) 
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The double ratio Rwa/ Rdq docs not depend on decay 
constants, and in our model we obtain Rd*o/ Roa = 
0.7955. The experimental value is given by Rd*q/Rdo = 
1.50 ± 0.75. Our result is small compared to the central 
experimental value but we are compatible within la. In 



the case of the meson Z?si(2460) we obtain 



Rdi = 0.7040 X 
Rd'i = 1.0039 X 



/Dsi(2460) 
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and for the double ratio Rd-'i/Rdi we get 1.4260, which 
agrees weh with the experimental result R^ti/ Rjji = 
1.32 ± 0.43. Finally, for the meson Dsi(2536) we obtain 



Rdi' = 0.6370 X 
Rd'I' = 0.9923 X 
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and for the double ratio Ro'V /Rdv, our value is 1.5578 
which in this case is 2cr above the experimental one, 
0.90 ±0.27. 

The quality of the experimental results does not al- 
low to be very conclusive as to the goodness of fac- 
torization approximation. But one can conclude from 
Eqs. (23), (24) and (25) that one cannot ignore, as done 
when using the infinite heavy quark mass limit, phase 
space and weak matrix element corrections. 

The decay constants of pseudoscalar and vector mesons 
in charmed and charmed-strange sectors are given in Ta- 
ble III. We compare our results with the experimental 
data and those predicted by different approaches and 
collected in Refs. [8, 34]. Our original values are those 
with the symbol (f). The decay constants of vector 
mesons agree with other approaches. In the case of the 
pseudoscalar mesons, the decay constants are simply too 
large. The reason for that is the following: Our CQM 
presents an OGE potential which has a spin-spin con- 
tact hyperfine interaction that is proportional to a Dirac 
delta function, conveniently regularized, at the origin. 
The corresponding regularization parameter was fitted 
to determine the hyperfine splittings between the n^So 
and rv^Si states in the different flavor sectors, achieving a 
good agreement in all of them. While most of the phys- 
ical observables are insensitive to the regularization of 
this delta term, those related with annihilation processes 
are affected because these processes are driven by short 
range operators [35]. The effect is very small in the '^5'i 
channel as the delta term is repulsive in this case. It is 
negligible for higher partial waves due to the shielding by 
the centrifugal barrier. However, it is sizable in the ^5*0 
channel for which the delta term is attractive. 

One should expect the wave functions of the 1^5*0 and 
I'^S'i states to be very similar [36] except for the very 
short range. In fact, they are equal if the Dirac delta 
term is ignored. The values with the symbol (|) in Ta- 
ble III are referred to the pseudoscalar decay constants 
which have been calculated using the wave function of 
the corresponding '^Si state. We recover the agreement 
with experiment and also with the predictions of different 



theoretical approaches. The foslfo and fD-Jfo* ratios 
are also shown in the last column of Table III. They are 
not very sensitive to the delta term and our values agree 
nicely with experiment and the values obtained in other 
approaches. 

Table IV summarizes the remaining decay constants 
needed for the calculation we are interested in. There, 
we show the results from the constituent quark model in 
which the 1-loop QCD corrections to the OGE potential 
and the presence of non-gg degrees of freedom in J^ = 1"*" 
charmed-strange meson sector are included. 

Refs. [37] and [38] calculate the lower bounds of the 
decay constants of I?*q(2317) and Z?si(2460) analyzing 

experimental data related to i? — > DD*j . Ref. [37] pro- 
vides the following lower limits: 



l«l|.fD:o(2317) 



58 - 83 MeV from B' decays 
63 - 86 MeV from 5° decays 



_/ 188^54 MeV from S- decays 
|ai| /d.,(2460) - < ^52t^^ MeV from S" decays 



(26) 



and the authors of Ref. 



get 



|ai|/. 



Z3;o(2317) 



74 ±11 



|ai|/D,i(2460) = 166 ±20 



(27) 



where the parameter jai | ^ 1. Our results arc compatible 
with these lower limits. 

Our results for the decay constants clearly deviate from 
the ones obtained in the infinite heavy quark mass limit. 
In that limit one gets /d*(,(2317) = Id,, /d,i(2460) = /r>* 
and /d^i(2536) = 0, results that lead to a strong dis- 
agreement with experiment for the decay width ratios in 
Eqs. (1) and (2). That was already noticed in Ref. [15], 
where the authors, using the experimental ratios, esti- 
mated that /D*g(2317) ~ \fDs and /d*o(2317) ^ /d,i(2460) 

instead. We obtain f Di„(2d.n) I f d , = 0.36, /d*^(2317) ^ 
0-72/d^i(2460) and /d,i(2536) = 59.176 MeV, the latter 
being small compared to the others but certainly differ- 
ent from zero. 

Finally, we show in Table V our results for the ratios 
written in Eqs. (1) and (2). The symbol (*) indicates that 
the ratios have been calculated using the experimental 
pseudoscalar decay constant in Table III. We get results 
close to or within the experimental error bars for the 
£'*o(2317) meson, which to us is an indication that this 
meson could be a canonical cs state. The incorporation 
of the non-qq degrees of freedom in the J^ = 1"*" chan- 
nel, enhances the jq = 3/2 component of the £'si(2536) 
meson and it gives rise to ratios in better agreement with 
experiment. Note that this state is still an almost pure 
qq state in our description [26]. 

The situation is more complicated for the Dsi(2460) 
meson. The probability distributions of its ^Pi and ^Pi 
components are corrected by the inclusion of non-gg de- 
grees of freedom, the latter making a ~ 25% of the wave 
function. In our calculation, only the pure qq component 



Approach 


fn (MeV) 


Id. (MeV) 


foJfD 


Ours 


297.019'+^ 


416.827^' 


1.40't' 




214.613") 


286.382**' 


1.33«' 


Experiment 


206.7 ±8.9 


257.5 ±6.1 


1.25 ±0.06 


Lattice (HPQCD+UKQCD) 


208 ±4 


241 ±3 


1.162 ±0.009 


Lattice (FNAL+MILC+HPQCD) 


217 ±10 


260 ± 10 


1.20 ±0.02 


PQL 


197 ±9 


244 ±8 


1.24 ±0.03 


QL (QCDSF) 


206 ± 6 ± 3 ± 22 


; 220±6±5± 11 


1.07 ±0.02 ±0.02 


QL (Taiwan) 


235 ± 8 ± 14 


266 ±10 ±18 


1.13 ±0.03 ±0.05 


QL (UKQCD) 


210±10l}^ 


236 ± 8lil 


1.13±0.02lH2 


QL 


211 ± 141^2 


231 ± 12tl 


1.10 ±0.02 


QCD Sum Rules 


177 ±21 


205 ± 22 


1.16 ±0.01 ±0.03 


QCD Sum Rules 


203 ± 20 


235 ± 24 


1.15 ±0.04 


Field Correlators 


210 ±10 


260 ± 10 


1.24 ±0.03 


Light Front 


206 


268.3 ±19.1 


1.30 ±0.04 


Approach 


fo' (MeV) 


fm (MeV) 


foi/fD' 


Ours 


247.865^' 


329.441*+' 


1.33(t' 


RBS 


340 ± 22 


375 ± 24 


1.10 ±0.06 


RQM 


315 


335 


1.06 


QL (Italy) 


234 


254 


1.04 ±0.01^4 


QL (UKQCD) 


245 ± 20I2 


272 ± 161^0 


1.11 ±0.03 


BS 


237 


242 


1.02 


RM 


262 ± 10 


298 ± 11 


1.14 ±0.09 



TABLE III. Theoretical predictions of decay constants for pseudoscalar and vector charmed mesons. The data have been taken 
from Ref. [8] for pseudoscalar mesons and from Ref. [34] for vector mesons. PQL=Partially-Quenched Lattice calculation, 
QL=Quenched Lattice calculations, RBS=Relativistic Bethe-Salpeter, RQM=Relativistic Quark Model, BS=Bethe-Salpeter 
Method and RM=Relativistic Mock meson model. 



Meson fo (MeV) y/M^fo (GeV"^/^) 



0:0 (2317) 
i:>si(2460) 
Dsi(2536) 



118.706 
165.097 
59.176 



0.181 
0.259 
0.094 



TABLE IV. Decay constants calculated within the CQM in- 
cluding 1-loop QCD corrections to the OGE potential and 
non-gg structure in channel l""". 



of the Dsi(2460) meson has been used to evaluate the 
T{B -^ D(*)Z?^i(2460)) decay width. The values we get 
for the corresponding ratios in Eqs. (1) arc lower than 
experiment. 



V. CONCLUSIONS 

We have performed a calculation of the ratios shown 
in Eqs. (1) and (2) working within the framework of the 
constituent quark model and in the factorization approx- 
imation. These ratios have been recently reported by 
the Belle Collaboration in Ref. [7] and provide impor- 
tant information to check the structure of the £'*o(2317), 
1)^1(2460) and i:>si(2536) mesons. 

The strong disagreement found between the heavy 
quark limit predictions and the experimental data moti- 



vates the introduction of the finite c-quark mass effects, 
which is done easily in the context of the constituent 
quark model. 

In the heavy quark limit and in the factorization ap- 
proximation, the ratios should be of order one for the 
-D|o(2317) and -Dsi(2460) mesons and very small for 
L'si(2536). In contrast, the experimental data indicate 
that the decay pattern of the Dsi(2536) follows the ex- 
pectations, but that is not the case for the Z?*q(2317) and 
£>.i(2460). 

The mass of the £'*o(2317) meson is lowered towards 
the experimental value with the inclusion of the 1-loop 
corrections to the OGE potential. We also obtain ratios 
compatible with the experimental data. Our results in- 
dicate that this meson could be described as a canonical 
cs state. 

We incorporate the non-gg degrees of freedom in the 
J^ = 1+ channel. The i3si(2536) meson remains almost 
a pure qq state and its jq = 3/2 component is enhanced. 
This gives the correct ratios for the Dgi (2536) meson. To- 
gether with the properties calculated in Refs. [26, 27], this 
is to us evidence of a good description of the Dsi(2536) 
meson. 

The Dsi{2A60) has a sizable non-qq component which 
contributes to the decays under study. This contribu- 
tion has not been calculated, as goes beyond the scope 
of this work. We have computed the ratios considering 
only the contribution coming from the qq structure of 





X = Dlo (2317) X = D,^ (2460) 
The. Exp. The. Exp. 


X = Dsi{2^m) 
The. Exp. 


B{B -^ DX)/B{B -^ DDs) 


0.19(*> 0.10 ±0.03 


- 


B{B ^ D''X)/B{B -^ D'Ds) 


O.IS'*' 0.15 ±0.06 


- 


B{B -^ DX)/B{B -^ DD;) 


0.177 0.44 ±0.11 


0.021 0.049 ±0.010 


B{B -> D*X)/B{B -^ 0*0*) 


0.252 0.58 ±0.12 


0.032 0.044 ±0.010 



TABLE V. Ratios of branching fractions for nonleptonic decays B — >■ D^^' D*J . 



the I?si(2460) meson. The ratios arc a factor 2 below 
the experimental ones. 
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Appendix A: DECAY CONSTANTS 

In our model, and due to the normalization of our nonrelativistic meson states, the decay constants are given by 



/m(o-) = 
/m(o+) = 
/a/(i-) = 
./'m(i+) = 



niM(Q-) 



mM(o^ 



ihf- 



mM{i- 



mM(n 



(0|Ji^^^(0)|M(0-),0) 



0|j4^^(0)|Af(0+),0) 



(0|4f (0)|Af(l-),0 0) 



(0|jif (0)|M(l+),00), 



(Al) 



for pseudoscalar, scalar, vector and axial-vector mesons, respectively. JyP^ and J^^^ are the vector and axial-vector 
charged weak current operators for a specific pair of quark flavors /i and /2. 

The corresponding matrix elements are given by 



(0|j{f (0)|A.f(0-),0) =^^ fd\p\ b1^0(X"(b1)t/- ^'''^''"'^^^'''^^'^"'"'^^ 



AEf,{-p)Ef,{p) 



\P\ 



(% (-P ) + n^/i )(% (P ) + ^h 
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.V3 



1 1 



l3liMiO+)), 



' JEf, i-p) + mf, ) {Ef, jp) + mf, ) 

4ii;/,(-p)i?;,(p) 



Ef^{p)+mf^ Ef,{-p) 



*/i 



V3 



(o|4Vno)|A/(i-,i = o),oo) = -^ / d\p\\p\'^f^'X 



X 1 + 



2i(M(l-,L=0))(|.|)^ 



l {Ef,{~p)+mf,){Ef,ip)+mf,) 
AEf,{^p)Ef,{p) 



3iEfA-p) + mf,){Ef,{p) + mj,) 



(0|4f (0)|M(1-,L = 2,0)0) =- 



2 /2 
^V3 



d\p\ ^ffjl-'^='^\\p\)Mf^^P^ + ^h)iEfA-p) + rufj 



AEi,{-p)Ef,{p) 



\P\ 



(-S/2 ip) + ™/2 ) (-E^/l (--P ) + "^/l ) ' 



7/2/1 



(0|4^3^^(0)|M(l+,5 = 0),00)=- / dblbTC/ 

1 



3;;(Af(l+,s=o)).|^|^ / (^/i (-P ) + mf, ){Ef,{p)+ ruf, ) 



4i?;,(-p)i?/,(p) 



. Ef^ i-p) + m/, Ef^ (p ) + m/. 



;o|.if^(o)|M(l^..l),oo)^-f/.b1b13^S:^-))(b1)J(^^^ 



(A2) 



Efi i~P) + ™/l Ef^ (p) + 771/, 

where Ef-^ and E'/, are the relativistic energy of the quarks. For 0^ and 0+ we have only one possible contribution. In 
the case of the J^ = 1 ^ meson we have two contributions coming from the two possible values of the relative angular 
momentum. For J = 1+ states and since C-parity is not well defined in charmed and charmed-strange mesons the 
wave function is a mixture of -^Pi and '^Pi partial waves and thus there are also two contributions. 



1. Decay constants in the heavy quark limit 

For the low-lying positive parity excitations, any quark model predicts four states that in the "^^^^Lj basis cor- 
respond to ^Pi, 3Po, ^Pi and ^P2. As charge conjugation is not well defined in the heavy-light sector, ^Pi and ^Pi 
states can mix under the interaction. 

In the infinite heavy quark mass limit, heavy quark symmetry (HQS) predicts two degenerate P-wave meson dou- 
blets, labeled by jq = 1/2 with J^ = 0+, 1+ (|l/2, 0+), |l/2, 1+)) and jq ^ 3/2 with J^ = 1+, 2+ (|3/2, 1+), |3/2, 2+)). 
In this limit, the meson properties are governed by the dynamics of the light quark, which is characterized by its total 
angular momentum jq = Sq + L, where Sq is the light quark spin and L the orbital angular momentum. The total 
angular momentum of the meson J is obtained coupling jq to the heavy quark spin, sq. 

A change of basis allows to express the above states in terms of the ^ ~^^Lj basis, by recoupling angular momenta, 
as 



(A3) 



|l/2,0+) 


= +\'Po) , 




11/2,1+) 


-+f/P^)U 


Jim 


13/2,1+) 


-^/P^)^^ 


flm 


13/2,2+) 


= +\'P2), 





where in the ^^'^^Lj wave functions wc couple heavy and light quark spins, in this order, to total spin S. 



11 





Do* (2400) 


Di(2420) 


D 


i(2430) 


Da* (2460) 


^Po 


+, 1.0000 


- 




- 


- 


'Pi 


- 


-, 0.5903 


~. 


0.4097 


- 


'Pi 


- 


+, 0.4097 


— . 


0.5903 


- 


'P2 


- 


- 




- 


+, 0.99993 


1/2,0+ 


+, 1.0000 


- 




- 


- 


1/2,1+ 


- 


+, 0.0063 


— . 


0.9937 


- 


3/2, 1+ 


- 


+, 0.9937 


+ : 


0.0063 


- 


3/2,2+ 


- 


- 




- 


+, 0.99993 




^:o(2317) 


D,i(2536) 


D, 


:i(2460) 


D:2(2573) 




Dio 


Dsi 




D'^i 


d:2 


^Po 


+, 1.0000 


_ 




_ 


_ 


'Pi 


- 


-, 0.7210 


— . 


0.1880 


- 


'Pi 


- 


+, 0.2770 


— . 


0.5570 


- 


'P2 


- 


- 




- 


+, 0.99991 


1/2,0+ 


+, 1.0000 


- 




- 


- 


1/2,1+ 


- 


-, 0.0038 


— . 


0.7390 


- 


3/2, 1+ 


- 


+, 0.9942 


— . 


0.0060 


- 


3/2,2+ 


- 


- 




- 


+, 0.99991 



TABLE VI. Probability distributions and their relative phases for the four states predicted by CQM in the two basis described 
in the text. In the 1+ strange sector the effects of non-qq components are included, see text for details. 



In the actual calculation the ideal mixing in Eq. (A3) between 'Pi and ^Pi states changes due to finite charm quark 
mass effects. Our CQM model predicts the mixed states shown in Table VI, which are very similar to the HQS states. 
This is expected since the c- quark is much heavier [rric = 1763 MeV) than the light (m„ = 313 Me V) or strange 
{nis = 555 MeV) quarks. Note that we also have mixing, even if small, between the ^P2 and '^^2 partial waves in 2+ 
mesons. This is due to the OGE tensor term. 

In Ref. [26] we have studied the J^ = 1+ charmed-strange mesons, finding that the J^ = 1+ Z3si(2460) has an 
important non-qq contribution whereas the J^ = 1+ Dsi{2536) is almost a pure qq state. The presence of non-qq 
degrees of freedom in the J^ = 1+ charmed-strange meson sector enhances the jq = 3/2 component of the Dsi (2536). 
This wave function explains most of the experimental data, as shown in Refs. [26, 27]. For this sector only the qq 
probabilities are given in Table VI. 

Here, we want to recover the results of the decay constants predicted by HQS. ^From Eqs. (A2), we only need to 
do the limit m/j — > oo 



Decay constants of Dg and D* in the heavy quark limit: 







(A4) 



EfAp) 



The Ds and _D* mesons belong to the doublet (|l/2, 0~), [1/2, 1^)) predicted by HQS. In such a limit the mesons 
are degenerate and the radial wave functions of the ^5*0 and '^Si components are the same. Therefore, we obtain 
fu^ = fo" in that limit. 



12 



Decay constants of D*q(2317) and -Dsi(2460) in the heavy quark hmit 



/r>*o(2317) 



1 



'^ V "^^;o(2317) 



d\p\\m'^jr'''''m)ji 



rrif^ 



Ef2 (P) Ef^ (p) + ruf^ 



/d»i(2460) 



-^ + 



"^I5si(2460) 



(A5) 



1 1 

6 TT 

2 1 



dwrn^'^jr^'^'-'m^l^ 



1 



^ll^|3l(-D=i(2460),='Pi) 



d\p\ \p 



^h,h 



(b1)Ji 



"^/2 

Ef2 (P) Ef^ (p) + m/, 
~~rn]~ 1 



Ef2 (P) Ef,_ [p) + m/. 



The L>*(,(2317) and D5i(2460) mesons belong to the doublet (|l/2, 0+), |l/2, 1+)) predicted by HQS. Again, in 
such a limit, the mesons are degenerate and the radial wave functions of the ^Pq, ^Pi and ^Pi components would 
be equal. Therefore, we would obtain /^. (2317) = /z)ai(2460) ™^ that limit. 

Decay constant of Z?si(2536) in the HQS: 



/m(1+, Dal (2536)) 



-^ + 



™I5,i(2536) 



^\\l I mmf.:}r'''''\mj^ 



3 77 

3^ 



771/2 



1 



rf|p I iP rc>: 



37(ri,i(2536),='Pi) 



(IpI)Ji 



£^/2 (P ) £^/2 (P ) + "^/2 
1 



(A6) 



m/. 



^/2 (P ) ^/2 (P ) + "I/2 



Therefore, considering that the radial wave functions of the ^Pi and '^Pi components are the same, /d^i(2536) is 
equal zero in the heavy quark limit. 
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